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We consider a single-species diffusion-limited annihilation reaction 2A — > with reactants con- 
fined to a two-dimensional surface with one arbitrarily large dimension and the other comparable 
in size to interparticle distances. This situation could describe reactants which undergo both longi- 
tudinal and transverse diffusion on long filamentous molecules (such as microtubules), or molecules 
that undergo truly one-dimensional translational diffusion (e.g. a transcription factor on DNA) 
but simultaneously exhibit diffusive behavior in a second dimension corresponding to a rotational 
or conformational degree of freedom. We combine simple analytical arguments and Monte Carlo 
simulations to show that the reaction rate law exhibits a crossover from one-dimensional to two- 
dimensional diffusion as a function of particle concentration and the size of the smaller dimension. 
In the case of a reversible binding reaction, the diffusion limited reaction rate is given by the Smolu- 
chowski expression, but the crossover is revealed in the statistics of particle collision histories. The 
results can also be applied to a particle-antiparticle annihilation reaction A + B — > 0. 



I. INTRODUCTION 

This paper studies a diffusion-limited single-species an- 
nihilation reaction 2A — > on a finite two-dimensional 
manifold for which one dimension is macroscopic in size 
and the size of the second is comparable to interparti- 
cle distances. Collisions between particles confined to 
this manifold lead to the formation of reaction products 
which are either inert or irreversibly dissociate from the 
manifold. Diffusion-limited annihilation reactions play 
important roles in a variety of physical and chemical sys- 
tems (e.g. pQ) and the kinetics of particle annihilation 
in low dimensions has been studied extensively in the 
context of particle-antiparticle annihilation (e.g. [2j [3]). 
It has been shown that correlations among particle po- 
sitions render the traditional Smoluchowski expression 
for the diffusion-limited rate incorrect for systems with 
fewer than three dimensions, and different rate laws have 
been established for one or two macroscopic dimensions. 
However, the effect of boundaries on low-dimensional re- 
action rates has not been thoroughly investigated, and it 
remains unclear what rate law applies for systems with 
extreme aspect ratios. 

While our results can be applied to particle- 
antiparticle annihilation, we were specifically motivated 
to consider reactions on manifolds with disparate di- 
mensions because many reactions in biological systems 
take place on one- or two- dimensional scaffolds, such as 
lipid-bilayer membranes, microtubules, or nucleic acids. 
It is becoming apparent that biological systems exploit 
low-dimensional manifolds to increase the rate at which 
molecules collide beyond the diffusion limited rate for 
molecules in three dimensions (3D). For example, exten- 
sive theoretical work (e.g. [3HE]) and in vitro experiments 
(e.g. [7HTT] and reviewed recently in [H]), and even in 
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vivo measurements |13[ I14j on transcription factors show 
that augmenting 3D diffusion with one-dimensional (ID) 
sliding of transcription factors on DNA molecules en- 
hances the rate at which they find specific target regions 
on DNA. Similar increases in reaction rates have been 
identified for membrane-bound G-protein signaling (re- 
viewed in |15j). and it has been proposed that the sliding 
of viral capsid protein molecules along an RNA molecule 
enhances the rate of capsid assembly [T&HT8] . 

We thus consider a physical system in which parti- 
cles diffuse on a filamentous molecule such as a long mi- 
crotubule or nucleic acid, and undergo bimolecular reac- 
tions to form complexes that either dissociate from the 
filament or become otherwise inert. Our analysis can 
also describe molecules that diffuse and react with spe- 
cific sites on the filament, such as a transcription fac- 
tor searching for a specific region of sequence. In the 
systems under consideration, the macroscopic dimension 
corresponds to the longitudinal direction along the fila- 
ment. The short dimension can correspond to the lateral 
direction if molecules undergo transverse diffusion (e.g. 
on a microtubule). Alternatively, for molecules that un- 
dergo true one-dimensional translational diffusion (e.g. 
along a DNA molecule), the short dimension can cor- 
respond to a rotational or conformational degree of free- 
dom whose value determines the reaction probability. For 
example, cross-sections for protein association reactions 
depend critically on the relative orientations of the col- 
liding reactants, and recent measurements indicate that 
a transcription factor has a relatively small probability 
of binding to its operator region each time it slides over 
it [T3]. We show that the reaction rate law in such a 
system exhibits a density-dependent crossover from one- 
dimensional to two-dimensional behavior as the size of 
the bounded dimension increases. 

The important implication of this result is that exper- 
imentally measured rate constants in such a system will 
depend on particle concentration for some system pa- 
rameters, with a form that itself varies with concentra- 
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tion. Our prediction of a crossover between one- and two- 
dimensional behavior also applies to particle-antiparticle 
annihilation reactions. Furthermore we show that for 
reversible diffusion-limited reactions, the particle corre- 
lations leading to alternative rate laws in low dimensions 
for irreversible reactions are reflected in the histories of 
particle-pair identities. 



II. THEORY 

A. Unbounded dimensions 

We consider a system of N particles which have dif- 
fusion constant D and annihilate upon collision; i.e. 
2 A — > 0. For the moment consider a system in d di- 
mensions, each of which has macroscopic size L. 

Traditionally the diffusion-limited rate 7Z is calculated 
under the "mean field" approximation that particle con- 
centration p — N/ L d is uniform in space and thus the bi- 
molecular collision rate is given by K = dp/dt = kp 2 /2. 
The rate constant k = WirRD with R the particle radius 
can be obtained from the Smoluchowski equation as de- 
scribed in the appendix. However, for fewer than three 
dimensions the inefficiency of the diffusion mechanism 
leads to amplifications of initial local density fluctuations 
and the reaction rate takes different, density-dependent 
forms. This was shown using dimensional analysis by 
Toussaint and Wilczec [3] and Krapivsky [2] as follows. 
In the mean field approximation it can be assumed that 
the rate constant depends on the particle radius and the 
diffusion constant, k — k(D,R), in which case dimen- 
sional analysis leads to k cx DR d ~ 2 with d the number 
of dimensions. To be physically realistic, the rate con- 
stant must increase with radius R, leading to the con- 
clusion that the mean field expression is only valid for 
d > 2. For d < 2 we must assume that the rate constant 
depends on the particle density rather than the radius, 
k = k(D,p), in which case dimensional analysis leads 
to k oc Dp~~ 1+2 / d and thus the reaction rate is given by 
TZ oc Dp 1+2 / d for d < 2. Since d — 2 is the marginal 
dimension, logarithmic corrections to the mean field ex- 
pression arise [2], which we derive approximately in the 
appendix. The dimension dependent reaction rates can 
then be summarized as 



ftiD = AiDp 3 
A 2 D P 2 
hx{B/p) 



K 2D = 



TZ 3B = 8nRDp 2 



(1) 
(2) 
(3) 



where A±, A 2l and B are unknown scaling constants, 
with Ai and A 2 determined by D and the collision cross- 
sections. 



B. Bounded dimension 

The expressions for the diffusion limited rates (Eqs. 
T][3) apply to a manifold that is unbounded in each di- 
mension. We now address the question of what happens 
on a manifold in which one of the dimensions is finite. 
We consider N particles in a 2D system with dimensions 
of size L and w with L >> w. Motivated by the physical 
systems described in the introduction, we consider the 
short dimension to be periodic, but the results are qual- 
itatively unchanged for reflecting boundary conditions. 
Without loss of generality we set the diffusion constants 
in each dimension equal to D; the size of the second di- 
mension can always be rescaled to achieve this situation. 

For sufficiently large w the particles are unaffected by 
the boundary conditions and we expect eq. |2]) to ap- 
ply, with p — N/Lw. However, below a threshold value 
of w, particles cover space in the short direction more 
rapidly than they collide. In this limit the short dimen- 
sion size will play no role, and the reaction rate will follow 
r = AiDpfjj with the effective one-dimensional concen- 
tration piD = N/L. We expect a crossover between the 
one- and two-dimensional rate laws at the threshold value 
of w. A mean field argument would suggest that the 
crossover occurs at the particle concentration for which 
the typical distance between particles is equal to the short 
dimension size, pio = w^ 1 ; however, the correlations be- 
tween particle positions described above change this esti- 
mate dramatically. Instead, noting that the frequency of 
collisions can never exceed that predicted by either equa- 
tion ([!]) or ([2]), the crossover can be determined as the 
short dimension size w and 2D particle concentration p c 
for which collision frequencies predicted by the ID and 
2D expressions are equal: 



Pc In 



B 



A 1 w 2 ' 



(4) 



This transcendental equation can be solved numerically 
for p given w, L, with the parameters B and A 2 /Ai de- 
termined independently from measured reaction rates. 



C. Reversible diffusion-limited reactions 

It is important to note that the equations described 
above apply to nonequilibrium irreversible reactions. In 
the case of a reversible diffusion limited reaction, where 
particles form dimer complexes with finite lifetimes, the 
mean-field prediction that reaction rate will be propor- 
tional to the square of the concentration must be correct, 
since the particle positions are distributed uniformly in 
the absence of energetic interactions to satisfy the Boltz- 
mann distribution [19) . However, the nature of diffusion 
in low dimensions and the resulting correlations described 
above are reflected in the history of particle collisions. 
Specifically, the rate of "novel collisions" , meaning colli- 
sions between pairs of particles which have never previ- 
ously interacted, will follow the forms given in @ or ((2|. 
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With particle tracking techniques and or single molecule 
fluorescence experiments [101 HU 031 E] it is possible to 
monitor the history of particle collisions in a reversible 
system. We show below that the expressions ([!]) and |2| 
are nearly recovered even if particle collisions are moni- 
tored only to first order, i.e., a collision between particles 
i and j is defined as "first-order novel" so long as parti- 
cle j is not the last particle that i was bound to and vice 
versa. 



III. SIMULATIONS 

To investigate the predicted rate laws ([T]) and ([2| and 
the density-dependent crossover size Q, we performed 
Monte Carlo simulations to study a diffusion-limited an- 
nihilation reaction on a two-dimensional rectangular lat- 
tice with periodic boundary conditions, with side lengths 
L and w measured in units of the lattice spacing a. In 
each Monte Carlo move a particle was chosen randomly 
and moved to one of its four neighboring lattice sites, 
chosen at random. If the move placed the chosen parti- 
cle on an already occupied lattice site, the two particles 
annihilated each other (i.e. they were deleted from the 
system). We measure distances in units of a and time in 
units of Monte Carlo sweeps, defined as the number of 
moves required for each particle to make on average one 
move in the long direction. Thus the diffusion constant 
was D = 1/2. 

Since the model was intended to investigate the effects 
of correlations and particle positions in a nonequilibrium 
system, a statistically appropriate initialization of parti- 
cle positions was essential. We began all simulations with 
100% of lattice sites occupied by particles, in which case 
there was only one possible arrangement and thus no a 
priori assumptions about the distribution were necessary. 
We then measured the annihilation rate as a function of 
the number of particles N(t) until the particle number 
decreased to near zero. 

We considered a second set of simulations to test the 
assertions corresponding to reversible collision-limited 
reaction rates. These simulations were performed as 
described for the annihilation simulations, except that 
when a Monte Carlo move would transfer a particle on 
to an already occupied lattice site, the initial positions of 
the two particles were swapped. This simulated forma- 
tion of a complex whose lifetime is short in comparison 
to the characteristic timescale for a particle to diffuse its 
diameter. The results pertaining to the form of the re- 
action rate law did not change when longer complex life- 
times were considered. In these simulations the previous 
reaction partner of each particle was stored in memory, 
and we recorded both the overall rate of collisions and 
the rate of first order novel collisions as defined above. 
These simulations were performed for a range of particle 
numbers N. 

In most simulations we set the long dimension size L = 
10 6 , but to enable sufficient statistics at low densities 
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FIG. 1. Simulations for all short dimension sizes w 
demonstrate a crossover from 2D to ID behavior, and 
that annihilation rates in the ID regime do not de- 
pend on w. Color online. Annihilation rates as a function 
of the effective one-dimensional concentration pm = N/L are 
shown for three representative values of w. The red solid lines 
and blue dashed line respectively show fits to the 2D (Eq. S) 
and ID (Eq. Q) rate laws. 



we performed additional simulations with L = 10 7 . We 
varied w in the range w e [5,500]. We performed 10 
- 100 simulations for each value of w. The simulations 
were optimized to efficiently handle large systems with 
concentrations ranging from 100% to near 0%. 



IV. RESULTS 

Simulations for all sizes w of the short dimension 
demonstrate a crossover from the 2D diffusion-limited 
rate law at high concentrations (i.e. large values of N) to 
ID scaling at low concentrations. Furthermore, the anni- 
hilation rates at different values of w below the crossover 
concentration collapse onto a single curve, consistent 
with the expectation of one-dimensional behavior. To il- 
lustrate these observations, annihilation rates as a func- 
tion of the effective one-dimensional concentration piD 
are shown for three short dimension sizes in Fig. [I] We 
see that the measured reaction rates at small or large 
densities agree very well with the theoretical rate laws 
for one and two dimensions, |l]) and (2|) (shown as solid 
red and dashed blue lines, respectively). The data was 
fit to the theoretical expressions independently at each 
value of w, with fit parameters Ai for TZiu arid A 2 for 
72.2D- The parameter B for 7?-2D is also unknown, but 
we held it constant at B = 1 to reduce the number of 
fit parameters. Further evidence that the theories match 
the simulation data is shown in Fig. [2j where we see that 
the independently fit values of A\ and Ai are constant 
within error over a wide range of w. 

Near the crossover, the measured rate is smaller than 
that predicted by either expression, as assumed in the 
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FIG. 2. Fit parameter values are independent of w. 

Color online. The fitted values for Ai and A2 are shown 
as functions of short-dimension width w. They are approxi- 
mately constant over the simulated range of w, except at the 
extremes where statistics are limited. The constant values 
for Ai and A2 used for the global fit in Fig. [3] are shown 
respectively as a solid red line and a dotted blue line. 
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FIG. 3. Measured crossover densities match theoret- 
ical predictions. Color online. The measured crossover 
concentration p c corresponds well to Eq. Q solved using con- 
stant values of A\ = 12.6 and A2 = 18.4 (solid red line, see 
Fig. [2| over the entire range of short-dimension sizes w. The 
simulated crossover densities (+ symbols) are obtained from 
averages over many independent simulations at each value of 
w as described in the text. If the logarithmic correction to 
2D behavior is neglected, the fit worsens (the best attempt is 
shown as a blue dotted line). 



theoretical analysis. The fact that the ID and 2D theo- 
retical rate laws fit the data for small and large w respec- 
tively, with essentially constant values of Ai and A 2 , then 
ensures that the theoretical prediction for the crossover 
density, p c in Eq. Q, will provide a good match for the 
data. To illustrate this, we define the simulated crossover 
densities as the densities for each value of w at which the 
rate laws fit respectively to small and large N intersect. 



FIG. 4. Reversible reaction rates follow mean-field 
scaling, but first-order novel collision rates approach 
the irreversible scaling. Color online. Collision rates are 
shown for a system in which particles swap places rather than 
annihilate upon colliding, with w — 10 over a range of effective 
one-dimensional concentrations pm well below the crossover 
density. The overall collision rate (red solid line) obeys the 
mean-field approximation for two dimensions (□). However, 
the rate of "first-order novel" collisions (0), defined in the 
text, nearly restores the ID irreversible collision rate scaling 
(dashed blue line). 



In Fig. [3] we then compare these simulated crossover den- 
sities to the the theoretical prediction Eq. Q calculated 
using constant values of A\ = 12.6 and A2 = 18.4. We see 
that the correspondence between the simulated and theo- 
retical crossover densities is excellent over two decades of 
short dimension sizes. We also show in the figure the best 
fit for the crossover densities if the logarithmic correction 
factor for the two-dimensional regime is not included in 
Eq. j2). This expression is slightly but consistently less 
accurate. 



Reversible reaction rates 



Reaction rates for the reversible model in which par- 
ticles swap rather than annihilate upon colliding were 
measured at w = 10 over a range of densities well be- 
low the crossover density p c = 2.0 x 10 -3 . In contrast 
to the irreversible annihilation model, which behaved ac- 
cording to the one-dimensional diffusion limited rate law 
at these densities, the overall reaction rates for this re- 
versible system are well fit by the mean field rate ex- 
pression 1Z = Dp 2 (Fig. [i]). Interestingly though, the 
one dimensional law is nearly recovered when the rate 
of first-order novel collisions (defined in section II C) is 
plotted as a function of density. 
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V. CONCLUSIONS AND OUTLOOK 



1. Fick's Law approach 



The results of our simulations demonstrate that 
diffusion-annihilation reactions on two-dimensional man- 
ifolds with large aspect ratios transition from one- 
dimensional behavior to two-dimensional behavior as the 
size of the shorter dimension crosses a threshold value. 
The simulation results are well fit by our theoretical ex- 
pressions, and confirm the significance of the logarith- 
mic correction that arises in the diffusion-annihilation 
rate law in two-dimensional systems. Furthermore, our 
simulations of the reversible reaction system show that 
particle correlations and the resulting transition to one- 
dimensional behavior at extreme aspect ratios are re- 
flected in the history of reaction partners. 

While we consider one macroscopic dimension and one 
small dimension, extension of our results to other situ- 
ations is straightforward. For example, in the case of 
one macroscopic dimension with length L and two small 
dimensions with sizes w and h (e.g. a particle diffusing 
within a narrow tube, or a particle undergoing transla- 
tional diffusion on a filament and two rotational dimen- 
sions) and a particle density p — N/Lwh, we expect a 
crossover from one-dimensional to three-dimensional be- 
havior at density p = R/d 2 wh. Similarly for diffusion and 
reaction of particles on a membrane, a single rotational 
dimension leads to a crossover between two-dimensional 
and three-dimensional (i.e. mean field) behavior. The 
important observation is that the crossover size of the 
"small" dimensions can become quite large for sufficiently 
small particle densities. 

Finally, extension of our results to particle-antiparticle 
annihilation, A + B — >• 0, is straightforward. While this 
reaction and the single species annihilation we consider 
here yield different kinetics for inhomogeneous starting 
conditions [2], the crossover relationship we predict is 
unchanged except within scaling constants. 
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Appendix A: Diffusion limited rates in one and two 
dimensions 



In this section we obtain scaling laws for the diffu- 
sion limited rate in one and two dimensions by two ap- 
proaches. 



First, we follow the classical solution for the diffusion- 
limited rate in three dimensions. In this approach, the 
rate constant for bimolecular collisions of particles each 
with radius R and diffusion constant D is mapped on to 
the problem of calculating the flux of point particles with 
diffusion constant D' — 2D to an infinitely absorbing 
sphere with radius R s — 2R [2U]. The concentration 
of point particles is held fixed at poo at a radius b > 
R. We then seek the steady-state flux of point particles 
to the sphere by solving the steady-state solution of the 
diffusion equation with p(r) the concentration of point 
particles at position r: 



V 2 p(r) =0 



(Al) 



with p(R s ) = 0, p(b) — poo and r — |r|. In three dimen 
sions with spherical symmetry this yields 

1 - R s /r 



p(r) 



l-Rs/b' 



(A2) 



The reaction rate constant k 3 D is then given by the area 
of the absorber, 4ttR 2 , times the inward flux, —J(R), to 
give 



ho = 4:irR s pooD 



1 



(A3) 



1 

R s /b' 

The usual expression for the diffusion limited rate is 
obtained by setting b — oo, substituting D' = 2D 
and R s = 2R and assuming the mean-field expression 
n 3D = k 3DPoo /2 = 8nRDp 2 . 

Following the same approach in one or two dimensions 
yields (assuming b^> R): 

AirDpoo 
ln(b/2R) 

--2D P oo/b 



km =: 



hi 



(A4) 



One immediately sees that the classical approach for 
three dimensions breaks down because there is no steady- 
state solution in the limit b — > oo. The situation can be 
rescued, however, by realizing that within the mean field 
assumption the bulk concentration must be established 
within a distance b = Poo from any molecule with d 
the number of dimensions. Substituting this relationship 
into ( |A4| ) yields 

8irDpoo 



«2D 



ln(l/4 Poo i? 2 ) 
fcio =2Dpl 



and thus 



n 



2D 



n 



ID 



-AirDpl 
--Dpi 



1 



ln(l/4 Poo i? 2 ) 



(A5) 

(A6) 
(A7) 



which are nearly equivalent to the expressions obtained 
by the more sophisticated quasi-static analysis that fol- 
lows. 



G 



2. Quasi-static approximation 

In this approach we begin by considering the rate ex- 
pression for the diffusion-annihilation reaction 2 A — > 0, 
with the density of species at position r and time t given 
by p(r, t) and mean density Poo{t), particle size R and dif- 
fusion constant D. Following Ref. [2] we start with the 
fundamental expression for the reaction rate as r ~ p/T 
with T the typical time interval before a given particle 
experiences a collision. Now we consider a tagged particle 
and work in a coordinate system with origin at the cen- 
ter of the tagged particle. Neglecting correlations among 
nearby particles with each other, their density satisfies 
the time-dependent diffusion equation 



dp(r, t + r) 



= 2L>V 2 .p(r,t + r) 



(A8) 



with the factor of 2 because the tagged particle diffuses 
as well. 

In the quasi-static approximation ^ [5T] it is assumed 
that density correlations over time r with respect to the 
tagged particle are limited to the distance it sweeps out 



in space by diffusion, \j2Dr. Thus we solve ( |A8[ ) with 
boundary conditions p(r = 2R) = for r > and p(r > 
\/2Dt, t) = poo to obtain 

, N 21n(r/2i?) , lM 
P(r,t + r), P (0,t) M ^ /2R l y (A9) 

The collision time T is then calculated as the time at 
which the flux across the circle of radius r = 2R equals 
unity: 



s.DRl Mr = 2R,t + T) dT = 1 
or 



(A10) 



Inserting (IA9J) into (A10) gives to leading order T 



ln(l/327T/Ooo i? 2 )/87rZ3p 00 . We then obtain the reaction 
rate 



TZ 2 u = 4ttZV 



1 



lnld^Trpi? 2 ) ' 



(All) 



The same analysis for d = 1 yields ( A7 ) 
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